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We provide a new perspective on the pseudogap physics for attractive fermions as described 
by the three-dimensional Hubbard model. The pseudogap in the single-particle spectral function, 
which occurs for temperatures above the critical temperature T c of the superfluid transition, is 
often interpreted in terms of preformed, uncondensed pairs. Here we show that the occurrence 
of pseudogap physics can be consistently understood in terms of local excitations which lead to a 
splitting of the quasiparticle peak for sufficiently large interaction. This effect becomes prominent at 
intermediate and high temperatures when the quantum mechanical hopping is incoherent. We clarify 
the existence of a conjectured temperature below which pseudogap physics is expected to occur. 
Our results are based on approximating the physics of the three-dimensional Hubbard model by 
dynamical mean field theory calculations and a momentum independent self-energy. Our predictions 
can be tested with ultracold atoms in optical lattices with currently available temperatures and 
spectroscopic techniques. 
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I. INTRODUCTION 

The analysis of peaks in the single particle spec¬ 
tral function, measured, for instance, by photoemis¬ 
sion experiments in solids or radio frequency (RF) spec¬ 
troscopy for ultracold atoms, provides important infor¬ 
mation about correlation effects in interacting quantum 
many-body systems. In the limit of weak interactions 
the spectral function displays peaks close to the ener¬ 
gies of the free fermion energy-momentum distribution 
and as such directly represents single-particle properties. 
At finite temperature and energies away from the Fermi 
surface, peaks are broadened and the width is indicative 
of interaction effects, which open up decay channels. If 
spontaneous symmetry breaking occurs below a certain 
temperature, such as in a superconductor below T c , the 
single-particle excitations become gapped out and shifted 
by an amount A, the superconducting gap. 

A more peculiar behavior is that of excitations be¬ 
ing gapped out (or suppressed) even though no obvious 
symmetry breaking and thermodynamic ordering transi¬ 
tion occurs. This is often referred to as pseudogap (PG) 
physics. The spectral gap can look very similar to a gap 
due to symmetry breaking at finite temperature. There¬ 
fore, it can be difficult to clarify the origin of PG physics 
and to distinguish whether it is due to some hidden or¬ 
der or a different effect. A very prominent example of 
such physics is provided by the experimental observa¬ 
tions in the hole doped copper-oxide high temperature 
superconductorswhere a relatively large part of the 
phase diagram is occupied by such a PG behavior. This 
phenomenon has attracted an enormous amount of atten¬ 
tion, however, there is currently no consensus about the 
physical origin of the this PG for the cuprates, and differ¬ 
ent scenarios have been invoked as an explanation. These 
include hidden order,® spin-fluctuations,® phase fluctu¬ 
ations and preformed pairsP 6 and the interplay with 


charge fluctuations. 7 

Here we focus on a conceptually simpler situation 
where PG physics has also been reported and that is 
for systems of fermions with locally attractive interac¬ 
tions. In situations without nesting the dominant insta¬ 
bility at low temperature is superconductivity and, cor¬ 
respondingly, pairing processes are expected to be most 
relevant. In particular, the crossover from weak coupling 
Bardeen, Cooper, and Schrieffer 8 ! (BCS) theory of super¬ 
conductivity to strong coupling Bose-Einstein condensa¬ 
tion (BEC) of pairs has been studied extensively and is 
a classical problem in condensed matter physics.® 12 In 
the last decade it has attracted renewed interest due to 
experimental realizations with ultracold fermions. Su¬ 
perfluidity has been reported in such systems,®^®® also 
in the case where the fermions are confined to an optical 
latticed Moreover, based on RF spectroscopy PG signa¬ 
tures have been reported for two 18 and three-dimensional 
systems without optical lattices 

Three non-exclusive concepts are usually invoked to 
discuss the origin of PG physics for attractive fermions: 
(i) Preformed pairs ; for intermediate coupling strength, 
pair formation without condensation is expected to oc¬ 
cur at a certain temperature T p which is larger than the 
superfhiid (SF) phase transition temperature T c . These 
preformed pairs can lead to PG formation as a certain 
binding energy is required to break the pair and re¬ 
solve a single fermion excitation. 12 20 21 1 This idea leads 
to a popular scenario for PG physics and is illustrated 
in a schematic phase diagram in Fig. [l] (ii) Pairing 
fluctuations above T c and their effect on single particle 
properties via a many-body self-energy can lead to PG 
physics.^ (iii) Phase fluctuations ; in a situation where 
fermions are paired one can imagine that a finite mag¬ 
nitude of the order parameter establishes locally, how¬ 
ever, no macroscopic coherent SF phase develops due to 
strong phase fluctuation.® In this situation the presence 
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FIG. 1: (Color online) Schematic phase diagram for the pre¬ 
formed pair scenario in the T — U plane with critical temper¬ 
ature T c , pairing temperature T p , Fermi liquid (FL) regime, 
and PG physics below T p (after Randeria 20 ). 



FIG. 2: (Color online) Phase diagram at half filling. We 
distinguish four different regimes: the superfluid phase (SF); 
a non-Fermi-liquid regime (NFL), which is separated into a 
region with PG in p{uj) and a region without PG (no PG); 
and a Fermi liquid regime (FL) below the temperature Tfl- 


of the ordering tendency related to a gap can then lead to 
PG signatures in the spectral functions. 2 ^ This behavior, 
which coincides with a small SF stiffness, is expected to 
be particularly pronounced in two-dimensional systems. 

Within one and the same calculation it is very difficult 
to obtain non-perturbative results and to include all rel¬ 
evant fluctuation effects. The purpose of this work is to 
contribute to a better understanding of the importance 
of particular effects in the lattice situation based on non- 
perturbative calculations. It is important to distinguish 
different setups when comparing the occurrence of PG 
physics for attractive fermions. First of all, dimensional¬ 
ity plays an important role in determining the strength of 
fluctuations, and in particular, the two-dimensional sit¬ 
uation has more pronounced fluctuation effects. More¬ 
over, results can differ in calculations for a model de¬ 
fined in the continuum and one on a lattice, such as the 
Hubbard model. A well known example is the T c curve 
which drops with the coupling strength on the lattice 
as 1/C7, whereas it approaches a constant in the contin¬ 
uum. Here we will analyze the attractive Hubbard model 
in three spatial dimensions. We will use the dynamical 
mean field theory (DMFT) approximation 231 to compute 
the self-energies and spectral function in the normal and 
SF phase. This approximation is non-perturbative in the 
interaction strength and therefore can describe very well 
the occurrence of preformed pairs. However, it does not 
include the effect of phase fluctuations (iii) and also does 
not include the effect of small momentum pairing fluctu¬ 
ations. The PG physics observed in our work can there¬ 
fore not be related to such effects. Phase fluctuations 
above T c are usually argued to be of minor importance 
for spectral properties in three dimensions. 

There is a substantial literature of previous work 
on BCS-BEC crossover and PG physics for attractive 
fermions, which however does not provide a clear and 


complete picture about PG physics. A popular approach 
is the diagrammatic T-matrix approximation)^* 24 * 25 * 
which captures well the effect of pairing fluctuations (ii). 
It was applied to the 2d Hubbard modePSI 27 and PG fea¬ 
tures have been found in the non-selfconsistent version, 21 
also in the continuum in twcP^ and three dimensions pUKD 
Selfconsistent T-matrix calculations for the 3d continuum 
model have found no PG in the spectrum. 32 33 However, 
in the two-dimensional case recently PG behavior was 
found There are also non-perturbative calculations, 
such as DMFT and quantum Monte Carlo (QMC) which 
found PG features in the continuum modeP^" 37 and for 
the Hubbard model at different filling factors. 38 44 The 
latter results were found to be in good agreement with a 
diagrammatic technique.^ It is worth noting that QMC 
techniques usually need to perform analytic continuation 
of imaginary axis data which can lead to uncertainties in 
results for spectral functions. DMFT studies, including 
cellular versions, for the attractive Hubbard model have 
been carried out in the normal phase)^®^ and in the 
broken symmetry phase l 42 l 43 l 5Q t El 
Our major results are the following: 

• For large enough coupling strength we find PG 
physics at temperatures T > T c . At half filling 
the PG remains for all temperatures above T c and 
therefore a pairing temperature T p (Fig. fTI) is not 
decisive to invoke the PG in the spectralfunction 
(see Fig. [2|. For different fillings the spectral func¬ 
tion is shifted due to the flattening of the Fermi 
function, such that the main suppression of spec¬ 
tral weight does not occur at uj = 0. 

• The occurrence of PG physics at high temperatures 
can be understood via split local excitations on lat¬ 
tice sites visible for strong enough interactions. 












3 


• PG physics in the spectral function is related to 
Non Fermi Liquid (NFL) properties of the self¬ 
energy (see Fig. |3j detailed definition below). 

• We demonstrate in detail how the PG transforms 
smoothly into the superconducting gap, when the 
temperature is lowered through T c (see Fig. |6|. 

The paper is organized as follows: In Sec. II we briefly 
describe our model and method. Sec. Ill discusses con¬ 
ceptual background about the occurrence of PG physics 
in relation to the self-energy. In Sec. IV and V we show 
results for spectra and self-energies at and away from half 
filling before concluding in Sec. VI. In the appendix we 
compare the DMFT-NRG calculations to iterated per¬ 
turbation theory and to T-matrix calculations. 


II. MODEL DEFINITION AND DMFT 
CALCULATIONS 

Our study is based on the three-dimensional attrac¬ 
tive Hubbard modelpBS] which in the grand canonical 
formalism reads 

H = + h.c.) - n^2n i(r - U y]n i)t n i4 , 

i,j,a ia i 

(i) 

with the chemical potential /i, the interaction strength 
U > 0 and the hopping parameters Uj. c\ a creates a 
fermion at site i with spin cr, and rq ?cr = cJ a Ci ?a . We 
take only a nearest neighbor hopping (—£), so that the 
non-interacting energy-dispersion relation in the three- 
dimensional cubic system is given as, k = (k x ,k y ,k z ), 

e k = -2t[cos(k x ) + cos (k y ) + cos(fc*)]. (2) 

The dispersion satisfies E- k = £ k - The corresponding 
density of states (DOS) is denoted by po(e). In the cal¬ 
culations we will use the hopping t and the bandwidth 
W = 12t as energy scales. 

The main method used to study the Hamiltonian 0 
is the dynamical mean field theory (DMFT) Within 
DMFT, we have to self-consistently solve a quantum im¬ 
purity model describing a single lattice site in the en¬ 
vironment of all other lattice sites. In order to calcu¬ 
late the self energy for this quantum impurity model, 
we mainly use the numerical renormalization group 
(NRG) which is able to calculate accurately expecta¬ 
tion values, Green’s functions, and self energies at zero 
and finit e tem perature d 57 * 581 also in the superconducting 
casef 53 ^ 59 ® Dynamical correlation functions are calcu¬ 
lated within the NRG by broadening of a large num¬ 
ber of discrete excitations in the Lehman representa¬ 
tion, and as such do not require analytic continuation. 
For our calculations, we choose a log-normal broaden¬ 
ing functiorP^^with unusually narrow and temperature- 
independent width, b = 0.3. One of the main reasons for 


this is to avoid a large transfer of spectral weight to high 
energies which can be particularly important at higher 
temperatures. Using this narrow broadening leads to ar¬ 
tificial oscillations in the spectra, which originate from 
the discretization of the bath in the NRG-calculation. 
In order to produce physical spectra, we finally smooth 
these oscillations by averaging over Ac o = 0.01 IF. This 
averaging is justified for the present purpose, because we 
do not expect very fine and sharp structures in our spec¬ 
tra on this energy scales to determine the physics of the 
PG. Furthermore, we carefully compared our NRG cal¬ 
culated spectra with iterated perturbation theory (see 
appendix). The latter technique does not require to 
broaden discrete excitations and provides therefore a use¬ 
ful test in a suitable parameter regime. 


III. FEATURES OF THE SPECTRAL 
FUNCTIONS AND SELF-ENERGY 

Before presenting the results of our calculations it is 
useful to discuss some basic features of the Green’s func¬ 
tions and self-energy, which will help us to better under¬ 
stand under which conditions PG physics occurs. In the 
literature PG physics is considered quite generally either 
for the integrated spectral function p{u) = ^ k p k (w), 

which is equivalent to the local spectrum pu(cj), or for 
fc-resolved spectra p k ( uj) close to the Fermi surface. We 
will consider both quantities in this paper. We note that 
a PG in one of these does not necessarily imply one in 
the other quantity. 

Let us first note that in the limit of high tempera¬ 
ture T W, U correlation lengths become small and the 
physics is dominated by local processes.^ This is seen, 
for instance, when we consider the bare single-particle 
propagator in imaginary time 

G%{r) = ~Jf ]Te ife ^e-^e^n F (&), (3) 

k 

where £, k = £ k — p and r G [0,/3). In the limit of high 
temperature, /3 = 1/T 0 and np(Cfc) —1/2. Then 

G?(r) in Eq. ^ becomes essentially local, ~ 5^, and 
spatial components with i ^ j vanish exponentially with 
length scale A ~ Quantum mechanical hopping is 

largely incoherent in this situation. This also means that 
DMFT based on a local self-consistent approximation can 
become very accurate in this high temperature limit. For 
instance, high temperature expansions for the three di¬ 
mensional Hubbard model agree well with DMFT cal¬ 
culations for thermodynamical quantities, and this re¬ 
mains to be t he ca se down to temperatures of the or¬ 
der T ~ pp/s 1^2164] Q ne should ? however, note that the 
self-energy of the three dimensional Hubbard model does 
not become completely k- independent even in the limit 
T^ooP 

What are the implications from this for the spectral 
function and the self-energy? The excitations in the 
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limit where local physics dominates are determined by 
the local part of the Hamiltonian, H\ oc = —fi i a n ia ~ 
V At half filling the chemical potential is 

fixed to (i = —U/ 2, and depending on the occupation 
n — 0,1,2 we have the energies E a = 0,17/2,0, respec¬ 
tively. Excitations in the spectral function have finite 
matrix elements for states where the particle number 
differs by one. Hence, in the spectral function excita¬ 
tion at energies A E = ±U /2 can be expected. The 
corresponding self-energy for the atomic problem reads, 
Yiii(uS) = 4 ’ w h ere T —)> 0. This implies ^-function 
peaks at ±{7/2 in the spectral function. In Sec. IV we 
will see that DMFT results at high temperature and large 
U are indeed of a similar form, Im£^(cj) = — 4 %) • 

Away from half filling the situation is more complicated, 
but similar features remain visible. If the peak in the 
self-energy is strong enough, we find in the spectral func¬ 
tion increased weight at uj = ±{7/2 and a suppression of 
spectral weight at the Fermi energy. These are the signa¬ 
tures of the PG in the integrated spectral function. For 
strong interactions this effect remains observable down 
to intermediate temperatures. In other words, the PG in 
p(uj) is related to the existence of Hubbard bands which 
are visible in the spectral function at all temperatures. 

We now discuss the appearance of a gap and PG in 
the momentum resolved spectral function p k (w). In the 
normal phase the Matsubara Green’s function reads 


Gfo (iuj n ) 


1 

iu) n 6s Ei{iujn) 


(4) 


where we have assumed a momentum independent self¬ 
energy as appropriate for DMFT calculations. The 
spectral function is obtained from analytic continuation, 
iuj n uj + ir/, 77 —> 0 , to yield 

Pk{U}) 7r[w-&-E*(w)] 2 + £'M 2 ' () 


We have separated real (R) and imaginary (I) parts of 
the self-energy. 

In the SF state we can include an explicit symmetry 
breaking term, Ag C , Ag C —» 0 for spontaneous symme¬ 
try breaking, and the non-interacting Green’s function 
matrix G k (iuj n ) has the form, 


G_l(iu n ) 1 = 


A 0 

1 - 


A° 


lUJ n 


6s 


(6) 


For the interacting system we introduce the matrix self¬ 
energy £ fc (icj n ) such that the inverse of the full Green’s 
function matrix G k (iuj n ) is given by the Dyson equation 


G k {iu n )~ l = G° k (iUn)- 1 - z k (iu n ). (7) 


The diagonal component of the fc-dependent Green’s 
function reads 


C2 ( 4 ^ 71 ) 

Cl,fc 6^n)C2,fc F12 (iu) n )Yj21 (iuJ n ) 


( 8 ) 




FIG. 3: (Color online) Schematic plot for real and imaginary 
part of the self-energy £(cj) in the FL (left) and NFL (right) 
regime. We also show the corresponding spectral function 
p(k F ,uj) which shows PG behavior in the NFL regime. The 
dashed diagonal line, uj, helps to identify the solutions of the 
equation uj — ReE(cj), and those positions roughly coincide 
with the PG peaks. 


with Cl ,k( iuj n) = u - 6c - Fn(icj n ), C2 ,k( iuj n) = iu n + 
Cfc — £22 (icj n ). The off-diagonal self-energy £ 12 ( 0 ;), in 
particular its real part, plays the role of a dynamic gap 
function, Re£i 2 (cj) ~ A. Therefore, low energy spectral 
excitations which correspond to uj = z(£ k — £^( 0 )) in 
the normal phase are shifted by the gap A to ±E k ~ 

±Z^J (Cfc - £^(0 )) 2 + A 2 , where z~ x = 1 - <? w £^(0) is 
the renormalization factor. Usually we associate the gap 
with a binding energy of pairs and hence we can interpret 
this energy shift as an energy required to break a pair and 
see a single-particle excitation. 

We now discuss the occurrence of a PG for momenta 
close to the Fermi surface in the situation where no off- 
diagonal self-energy is present. Thus consider k = k F 
(interacting Fermi surface) such that Cfc F — £^(0) = 0P^ 
Then we can write 


1 £ 7 (cj) 

7 r [w-S fl (o ;)] 2 + S J (a;) 2 ’ 


(9) 


where H R (cj) = ’E R (u>) — £ fl (0). Provided that S 7 (w) 

does not vary rapidly, we expect p kF (uj) to be peaked 

_ 

when the implicit equation uj = £ (uj) is satisfied. Ac¬ 
cording to our definitions there is always a solution to 
this equation for uj = 0. In a weakly interacting system 
at low temperature |£ 7 (u;)| usually has a local minimum 
at uj = 0 , 

Im£(u;) = —a(T) — buj 2 , (10) 


where a(T) —0 for T 


0. By 1 
ig- 0 ( 


0 and a, b > 0. By the Kramers- 
(left)]. Then 


Kronig relation cG£ (0) < 0 [see Fig 

_ 

the only solution of uj = £ (uj) is the one at uj = 0 . 
This is the Fermi liquid peak in the spectral function 
at uj = 0 with width ^ z|£ 7 ( 0 )| and weight z, where 


z~ x = 1 — cG£ (0). We define the low energy behavior 

in Eq. (10) as Fermi liquid (FL) regime. _ 

A PG is obtained with different behavior.If 
|£ J (u;)| possesses a local maximum at uj = 0 , 


—R, 


Im£(cj) = —a(T) + buj 2 , 


G k (uj) 


( 11 ) 
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_ 

then c^E (0) > 0. If the slope is large enough we will 

_ 

then encounter additional solutions of uj = E (uj) as can 
be easily seen graphically [see Fig. [3] (right)]. Whether 
this is the case depends on the interaction strength, filling 
fraction and temperature. Since |E 7 (cj)| is decreasing, 
we obtain a local minimum at uj = 0 in the spectral 
function and broadened peaks at finite energies. This 
means that the original peak at uj = 0 is split and hence 
we obtain a PG. Notice that a local maximum of |E 7 (c<;)| 
does not necessarily lead to a PG, if the self-energy is 
not large enough. In the following we call the low energy 
behavior of Eq. Non-Fermi liquid (NFL) behavior. 
As we have discussed above |E 7 (cj)| is typically maximal 
at u) = 0 at high temperature when the physics becomes 
dominated by local interactions. It is also directly visible 
in the phase space factor appearing in the second order 
perturbation theory in U (see appendix). Therefore, at 
high temperature we expect NFL behavior, and at low 
temperature we usually have FL behavior. We define 
the crossover scale as Tfl, i.e. where the behavior of 
E 7 (cj) changes from Eq. (11) to (10). In this picture PG 
behavior in p^ F (uj) occurs therefore as long as (i) U is 
large enough (^ W) and (ii) T > Tfl(E 7)- In particular, 
the PG is always present above T c if T c > Tfl(£7)- 


IV. PG PHYSICS AT HALF FILLING 


In this section we analyze results from the DMFT cal¬ 
culations for spectral functions and self-energies and fo¬ 
cus on the situation at half filling. An overview of the dif¬ 
ferent regimes as function of U and T is shown in Fig. [2] 

The phase diagram includes the SF phase and the 
regimes where the self-energy shows FL and NFL be¬ 
havior as defined in Eq. (10) and Eq. respectively. 
By performing calculations suppressing the SF phase be¬ 
low T c , we find that the boundary between FL and NFL 
regimes (not shown) is connected to the bipolaron transi¬ 
tion at T = 0, which is equivalent to the Mott transition 
for repulsive interactions. The NFL regime in the phase 
diagram is separated into a region for stronger interac¬ 
tions where we observe a PG in the integrated spectral 
function, and a region without PG (no PG) for weaker 
couplings. 

In the upper panels of Fig. [4] we show the interact¬ 
ing local DOS p( uS). At weak coupling and intermedi¬ 
ate temperatures, p(uj) very much resembles the non¬ 
interaction DOS, po(u), and the small self-energy does 
not have a pronounced effect. Although |ImE| is peaked 
at the Fermi energy for U = 0.4W at high temperatures, 
there is no PG structure in the DOS. In contrast, for 
larger interactions, U/W = 0.6, U/W = 1, we find at 
high temperatures a PG structure of two peaks at ±17/2 
and a suppression of the density of states at uj = 0. The 
behavior is more pronounced for larger interactions. In 
both cases the magnitude of the PG is clearly related 
to U. This structure is induced by the NFL peak in 
the |ImE| (lower panels). As discussed in the previous 


U/W=0.4 U/W=0.6 U/W=1.0 




FIG. 4: (Color online) Integrated DOS, pn(uj), and imagi¬ 
nary part of the self-energy, Sn(w), for different interaction 
strengths and temperatures. 


section III this result can be understood in terms of the 
local excitations dominating the physics at high tempera¬ 
ture. At weak and intermediate interaction strengths the 
system crosses over to a FL regime before T c is reached 
when decreasing the temperature. For V/W = 0.4 and 
U/W = 0.6, |ImE(u;)| exhibits a dip at the Fermi en¬ 
ergy at low enough temperature, T/W < 0.1, which is 
accompanied with a peak structure in the DOS. Such a 
change in the behavior of self-energy and DOS cannot 
be observed for strong coupling, where the PG structure 
exists for all temperatures above T c . At very low tem¬ 
peratures, the system is in the SF phase in all cases, 
which is characterized by a gap in the DOS, which co¬ 
incides with a dip in ImE(u;). So even though the two 
cases, U/W = 0.6, U/W = 1, in Fig. [4] look similar at 
high temperature (PG) and very low temperature (SF 
gap), they display a striking difference for intermediate 
temperatures. For the larger coupling strength the SF 
transition occurs from a PG state (see also Fig. [6|; in 
contrast for U/W = 0.6, the SF instability happens in 
the FL regime. 
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FIG. 5: (Color online) The local pair density for differ¬ 

ent temperatures and interaction strengths. The black arrow 
marks the transition from non-Fermi-liquid to Fermi-liquid 
behavior. The red arrow marks the transition into the SF 
phase. 


Further insights can be obtained by studying the be¬ 
havior of the double occupancy or local pair density, 
(rzqrzq), which is shown in Fig. [ 5 ] for different temper¬ 
atures and interaction strengths. 

Independent of the interaction strength, in the high 
temperature limit, T W, U, the pair densities ap¬ 
proach the non-interacting values n 2 , where n a is the 
density for one spin component, at half filling ( 71 ^ 74 ) = 
0.25. In the atomic limit, t = 0, the double occupancy 
can be easily calculated. At half filling, all atomic states 
are occupied with equal probability, so that the double 
occupancy reads 

= 2 + 2exp(-C7/(2T))' (12) 

At high temperature, T/W > 0.5, this formula agrees 
very well with the results in Fig. [5j demonstrating again 
that the physics at high temperature is dominated by 
local processes. 

Decreasing the temperature, {n^n^) increases due to 
the attractive interaction. This effect is stronger for 
stronger interaction. For interaction strengths U/W < 
0.8, we find a maximum before the system enters the 
SF phase at T c . This maximum appears to be corre¬ 
lated with the crossover temperature Tfl (black arrows) 
between FL and NFL behavior in the self-energy. The 
disappearance of the maximum in the pair density for 
interaction strengths U/W > 0.8 agrees with the van¬ 
ishing of the FL regime phase in the phase diagram. For 
U/W < 0.8, the pair density decreases when lowering the 
temperature below Tfl, but then increases again when 
entering the SF phase (arrow at T c ). For strong interac¬ 
tions ( U/W > 0.8) on the other hand, the pair density 
increases with decreasing temperature until T c is reached 
and then decreases. This agrees with the known fact 


that the superfluidity is driven by interaction energy gain 
for weak coupling, as opposed to kinetic energy gain for 
strongly coupled systems. 52 

With these insights we can comment on how our re¬ 
sults compare to the preformed pair scenario in Fig. [l] 
It is interesting to note that for very high temperatures 
the PG behavior in Fig. [4] does not change significantly 
anymore. In other words the PG persists and no T p for 
its appearance can be identified. This is the case even 
for temperatures where the pair density has decreased to 
values close to the non-interacting result. Furthermore 
we found PG behavior for the two cases U/W = 0.6, 
U/W = 1 at high temperature, but for intermediate tem¬ 
peratures ( T/W ~ 0.05) the case U/W = 0.6 shows FL 
behavior. In both cases we observe a strongly enhanced 
local pair density for such temperatures, which can be 
interpreted as a preformed pair state, however, the man¬ 
ifestation in the spectral function is different. Both of 
these observations are in clear contrast to the preformed 
pair scenario, where the existence of the PG behavior is 
linked to the presence of an enhanced pair density! 20 21 

In Fig. [6j we take a closer look at dynamic response 
functions close to the SF transition temperature T c . 

The plots in the upper part of figure show p(ui) for 
U/W = 0.6 and U/W = 1, which correspond to a tran¬ 
sition into the SF phase from the FL and NFL regime, 
respectively. The lower part of the figure displays the cor¬ 
responding diagonal and off-diagonal self-energies (real 
and imaginary parts). For the weaker coupling case, 
U/W = 0.6, at T > T c there is the usual FL dip in 
Im£n(u;) and the corresponding peak in p{uS). When 
the temperature is lowered through T c the off-diagonal 
self-energy becomes finite and a dip in p(uf) is induced. 
Very close below the transition temperature, the main 
effects for this come from ReE^cj). Lowering the tem¬ 
perature further, the amplitude of the diagonal part of 
the self-energies decreases without showing new features. 
As discussed in Sec. Ill the gapping out of excitation is 
dominated by contributions from Re£i 2 (cj). 

In the case of stronger interaction, U/W = 1, super¬ 
fluidity sets in the NFL regime with a PG at the Fermi 
energy. When lowering the temperature through T c , the 
off-diagonal self-energy becomes finite, but at first the di¬ 
agonal part of the self-energy remains nearly unchanged 
(The orange line, T/T c = 1, overlaps with the dark green 
line, T/T c = 1.1). On further reducing the temperature 
the off-diagonal self-energy increases substantially and 
Im£n(u;) is strongly reduced developing a FL dip at the 
Fermi energy. The gap in p(ui) changes smoothly from 
the PG with broad peaks separated by U to the sharper 
structures (coherence peaks) in the SF phase. It is in¬ 
teresting to note that the gap, if defined as the distance 
between the maxima is larger above T c in the PG regime 
than in the SF phase. One should also note that for low 
temperatures the gap becomes much more pronounced 
with a suppression of spectral weight at u) = 0 and as 
such is approaching a full gap in the limit T 0. 

A remarkable observation is that the qualitative behav- 
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U/W=0.6 U/W=1.0 



FIG. 6: (Color online) p(uS) and £(u;) close to the SF phase 
transition. We use the same legend for the self-energies as 
for the Greens functions. The transition temperatures are 
T c /W = 0.32 for U/W = 0.6 and T c /W = 0.35 for U/W = 
1 . 0 . 


ior of the off-diagonal part of the self-energy can change 
within the SF phase. Generally, |Re£i 2 (tcQ| approaches 
the mean field result U (c^c^) for \uj\ — ooP^ At weaker 
coupling (U/W = 0.6) and low temperature it is minimal 
for small c o. Decreasing the temperature, the anomalous 
expectation value increases and this is reflected in the 
results for £ 12 ( 0 ;). The cj-dependence can be understood 
at weak coupling from the effective interaction for induc¬ 
ing superfluidity, which possesses a repulsive component 
which is peaked for small However, when entering 
the SF phase from the PG regime at stronger coupling 
( U/W = 1), |Re£i 2 (cj)| first develops a strong maximum 
at uj = 0. When the temperature is lowered further this 
behavior continually reverts to the one of the weak cou¬ 
pling situation. The form of the spectral function changes 
at T c and there is a shift from the gap feature being in¬ 
duced by £n (u) (above T c ) to £ 12 ( 0 ;) (below T c ). The 
observed strong changes are related to this shift and a 
more thorough understanding requires further investiga¬ 
tion. 


We now turn our attention to features in the momen¬ 
tum resolved spectral function pk(w). A good overview of 
the behavior for different interactions and temperatures 
can be obtained in the intensity plots in Fig. [7J 

We show pfc(cj) for three interaction strengths [U/W = 
0.4 (upper panels), U/W = 0.6 (middle panels), and 
U/W = 1 (lower panels)] for three different tempera¬ 
tures [T/W = 0.2 (left), T/W = 0.08 (middle), and 
T/W = 0.01 (right)]. We also show the Fermi level 
(dashed line) and the non-interacting dispersion (full red 
line) as an orientation. At weak coupling, U/W = 0.4, 
the spectral function only displays a weak modification 
from the non-interacting result with certain broadening 
of the peaks and a minor shift of spectral weight. At 
low temperature the system is SF and excitations at k F 
are gapped out. Notice that the width of the Bogoliubov 
peaks at the gap edge is overestimated by our broadening 
procedure.^ 

For U/W = 0.6, we find similar features for pk(u) as 
what has been found for the integrated spectral function, 
p(cj), as far as the PG is concerned. At high tempera¬ 
tures we see a broadened dispersion similar but shifted 
from the non-interacting one. Spectral weight is sup¬ 
pressed at the Fermi energy such that PG features are 
realized at high temperatures. Curiously, this PG closes 
at intermediated temperatures ~ Tp l where the behav¬ 
ior of the self-energy changes. Below T c the spectrum is 
gapped again. Notice that band renormalization features 
appear somewhat weaker than at high temperatures. 

For U/W = 1 the NFL regime extends from high tem¬ 
peratures down to T c . The self-energy undergoes only 
very slight changes when decreasing the temperature in 
the NFL regime. Accordingly, the momentum-resolved 
spectral function for T/W = 0.2 and T/W = 0.08 (lower 
left panel and lower middle panel) are nearly the same. 
We observe a large PG around the Fermi energy; the 
spectral weight at the Fermi energy is very small. When 
entering the SF phase, gap features are visible and the 
dispersion changes in the vicinity of u) = 0. For this in¬ 
teraction strength, we observe a clear deviation between 
the non-interacting band structure and the interacting 
spectral function. In the SF phase we find a mirror or 
“shadow” band appearing as reflected from uj = 0. These 
bands can be understood due to a particle-hole doubling 
in the Nambu representation. This is an effect also ob¬ 
served in the antiferromagnetically ordered phase with 
zone doubling. 70 

In Fig. [8] we show particular cuts for p& F ( uj) as a func¬ 
tion of uj. Here we can see the PG features even more 
clearly. Similar as in the integrated spectrum, p{ uj), the 
PG is absent for the weak coupling case, U/W = 0.4, 
but present at high temperature for stronger interactions, 
U/W = 0.6 and U/W = 1. For U/W = 0.6 the PG disap¬ 
pears in the FL regime, whereas it remains for U/W = 1. 
We also show the real part of the diagonal self-energy. As 
discussed in Sec. Ill, the peak splitting in pk F {u) can be 
induced from non-trivial solutions of uo = Re£n(u;), and 
we have included a dashed line to see this graphically. 
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FIG. 7: (Color online) Momentum resolved spectral function for U/W = 0.4 (upper panel), U/W = 0.6 (middle panels) and 
U/W = 1 (lower panels). The temperature are T/W = 0.2,0.08,0.01 from left to right. The red line corresponds to the 
non-interacting system, the green line corresponds to the Fermi energy. 




-0.4-0.2 0 0.2 0.4 -0.6-0.3 0 0.3 0.6 -0.8-0.4 0 0.4 0.8 
co/W co/W co/W 


FIG. 8: (Color online) pk F {w) for k F m (7r/2, 7 t/2, 7t/2) for 
U/W = 0.4, U/W = 0.6 and U/W = 1 for different tempera¬ 
tures. 


As can be clearly seen, this condition is not satisfied in 
the weak coupling case. In contrast, at strong coupling, 
U/W = 1, the intersection points characterize the peak 


positions well. The spectral function changes in the SF 
phase (lowest temperatures), where the coherence peaks 
at the gap edge become visible. 


V. PG PHYSICS AWAY FROM HALF FILLING 

So far we have focused on the situation at half filling 
where the discussion is somewhat simplified due to the 
particle-hole symmetry. In this section we show results 
for different filling factors (n < 1) to see how the PG 
behavior is affected. This is important for comparison 
with experiments with ultracold atoms, where due to the 
trapping potential no homogeneous filling fraction can be 
expected. 

In Fig. [9] results for p(uj) and Im£(u;) analogous to 
the ones in Fig. [4] for the half-filled case are displayed 
for n = 0.5 over a wide temperature range. The lowest 
temperature corresponds to a gapped SF state. Looking 
at the self-energies in the lower panel, we can clearly see 
that the classification into FL and NFL regions is still 
applicable and |ImE(u;)| can either show a double peak 
with dip in the vicinity of uo = 0 (FL) or a strong single 
peak (NFL). It is useful here to distinguish temperatures 

































































9 



FIG. 9: (Color online) The DOS and imaginary part of the 
self-energy for U/W = 0.6 and U/W — 1 for different tem¬ 
peratures. The filling of the system is fixed to n — 0.5. 

T/W <0.2, where features are close to uj = 0, and higher 
temperatures, where the NFL peak in |ImE(u;)| moves 
systematically to higher energies. In contrast to the half 
filled situation the case U/W = 0.6 does not show a clear 
PG in p(uj). However, for U/W = 1 the PG is clearly 
visible. For lower temperatures the minimum in p(uS) is 
close to uj = 0 and for higher temperatures it moves to 
higher energies together with the NFL peak in |ImE(u;)|. 
Notice, however, that the minimum in p(u) and the peak 
in |ImE(u;)| do not coincide as they do for n = 1. The 
shift of the PG with temperature can be understood by 
recalling that at high temperature the Fermi distribution 
becomes flatter such that higher energies contribute to 
the particle number, n = f du)p(w)riF(w). To satisfy this 
relation at higher temperature the spectrum has to be 
shifted. 

The PG structure in p{uj) at elevated temperature can 
still be understood from the local picture. For n < 1, we 
can write p = —U/2 — Ap (U > 0) assuming Ap > 0, 
and the atomic energies are E a = 0,?7/2 + A/x,2A/x. The 
partition function reads 

Z = 1 + 2e~ /s{u / 2+Afi) + e-? 2 ^ (13) 



FIG. 10: (Color online) The DOS and imaginary part of the 
self-energy for U / W = 0.6, and U/W = 1 for different fillings. 
The temperature of the system is T/W = 0.05. All data 
shown corresponds to the normal phase. 

There are now excitations at = U/2 + A p and c = 
— U/2 + Ap with generally asymmetric weights 

w + = -I [1 + e -P( u / 2 + A ri ], (14) 

Z 

and 

W_ — p-Q-PW/ 2 +A/x)]^ (15) 

Z 

respectively. 

Without showing explicit results we note that the pair 
density {n^n±) displays a similar temperature depen¬ 
dence for n = 0.5 to what was shown in Fig. [5} increas¬ 
ing from n 2 a at large T to larger values (maximal n/2). 
Therefore, similarly to the half filled case PG behavior 
can coincide with an enhanced pair density for large in¬ 
teractions and T c < T. However, we also find cases, e.g. 
U/W = 0.6, T/W = 0.05, with enhanced pair density 
((n^ni) ~ 0.17) and no PG behavior, in contrast to the 
expected relation in the preformed pair scenario. 

In order to get an insight to overall trends, we compare 
several different fillings in Fig. [lOj We show p(w) and 
ImE(u;) for U/W = 0.6 and U/W = 1 for low tempera¬ 
ture, T/W = 0.05, in the normal phase. For U/W = 0.6 
p(uS) exhibits a FL dip in ImE(o;). It is clearly visible, 
even for n = 0.1, that the self-energy does not change 
its structure when reducing the filling further. The fre¬ 
quency dependence in this FL regime is relatively sym¬ 
metric with respect to uj = 0. The DOS, on the other 
hand, changes with n. While for n = 0.5 a clear peak 
close to the Fermi energy is visible in the DOS at low 
temperature, such a peak is hardly noticeable for n = 0.2, 
and it has disappeared for n = 0.1. The amplitude of the 
self-energy has become too weak to change the spectrum 
and we essentially see a shifted non-interacting DOS. 

For U/W = 1, we find a NFL peak in ImE(cj) for 
all fillings in Fig. [lOj We observe similar effects as for 
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FIG. 11: (Color online) Momentum resolved spectral function pk(w) for n = 0.5 from top to bottom U/W = 0.4, U/W = 0.6, 
and U/W = 1; and from left to right T/W = 0.2, T/W = 0.08, and T/W — 0.01. The red line corresponds to the non-interacting 
dispersion £&, the dashed green line corresponds to the Fermi energy. 


weaker interactions when reducing the filling as far as 
the strength of the self-energy is concerned. However, 
we find a clear PG structure in the DOS. Whilst the 
PG structure at low temperature is pinned to uj = 0, 
at high temperature the whole spectrum including the 
PG is shifted to high frequencies (see Fig. |9|. Note that 
at high temperature due to the flattening of the 

Fermi energy (lj = 0) does not play such an important 
role as it does for low temperatures. 

In summary, when analyzing p{ui) and ImE(cj) we can 
find similar features to the ones of the half filled situation 
and a PG appears for suitable parameters. However, de¬ 
pending on filling, temperature, and interaction strength, 
the occurrence of the PG may be limited. At high tem¬ 
peratures it can be shifted away from uj = 0, although 
it is still clearly visible in the spectrum. Moreover, the 
impact of the local Hubbard interaction becomes weaker 
for a system with smaller filling factor. 

Momentum resolved spectra for n = 0.5 and various 
values of T and U are displayed in Fig. EH Generally, 
the features are similar to the half filled case. For weak 
coupling ( U/W = 0.4) we find a shifted and broadened 
spectrum which shows a SF gap at low temperature. For 
intermediate coupling (U/W = 0.6) interaction effects 


are more visible in the spectrum, resulting in stronger 
band renormalization effects and shifts of spectral weight. 
However, in contrast to n = 1 no clear PG becomes vis¬ 
ible in ph(w). For U/W = 1, we see strong interaction 
effects and PG features at all temperatures above T c . We 
also clearly observe an asymmetry in the intensity, which 
is substantially lower for the uj < 0 part of the spectrum. 


Particular cuts along uj for momenta which satisfy 
£k F + ReE(0) = 0 are shown in Fig. 12 At weak cou¬ 


pling ( U/W = 0.4) the FL peak is gapped out when the 
temperature is lowered below T c . For intermediate cou¬ 
pling ( U/W = 0.6) above T c we find that the FL peak 
is shifted away from uj = 0 to higher energies. Also the 
coherence peaks below T c show some asymmetry due to 
self-energy effects. A clear PG is only visible for larger 
interactions, U/W = 1. The lower panel shows again the 
real part of the self-energy. In contrast to the situation 
at half filling, the peaks in pk F (w) are not well explained 
by the intersection, uj = ReH^). In this situation the 
variation of ImE(o;) is too strong invalidating the simple 
arguments of Sec. III. Nevertheless a NFL peak form of 
the self-energy is clearly important for the PG behavior. 
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U/W=0.4 U/W=0.6 U/W=1.0 



FIG. 12: (Color online) p kF (w) for U/W = 0.4, U/W = 0.6 
and U/W = 1 and different temperatures. Notice that the 
lowest temperature T/W = 0.01 is below T c in all cases. The 
lower panels show the real part of the self-energy including 
the line y — uj. 


VI. DISCUSSION AND CONCLUSIONS 


setup. 

Experiments with ultracold atoms in optical lattices 
provide an excellent platform to test our predictions. In¬ 
teractions can be tuned in a wide range by Feshbach res¬ 
onances, the lattices can be loaded with different filling 
factors and a temperature range T/W = 0.1 — 0.2 is rou¬ 
tinely accessible 16 . Integrated and momentum resolved 
spectra can be measured such that a direct comparison 
with our predictions is possible. Thus, we hope that our 
work will stimulate further efforts in this field which con¬ 
tribute to a better understanding of the intriguing PG 
physics. 
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We have analyzed the occurrence of PG features in 
the integrated and k -resolved spectral function of the 
three-dimensional attractive Hubbard model for different 
temperatures, interactions and filling factors. Properties 
of the spectral functions have been traced back to the 
characteristic behavior of the self-energy. We find PG 
behavior as long as the interaction U is large enough 
W) and the self-energy shows NFL behavior, i.e., 
T > Tfi,(U). Our results show marked deviations from 
the popular preformed pair scenario, where PG behavior 
is directly linked to the formation of pairs at a tempera¬ 
ture T p > T c : (i) We find that PG behavior persists up 
to large temperatures and is not bounded by some tem¬ 
perature scale T p . (ii) We find cases with a substantially 
enhanced pair density where no PG behavior occurs. The 
first effect is related to the fact that we are working with 
a lattice model, such that local excitations are always 
well defined and related to the chemical potential fi and 
U . This might be different in the continuum where it is 
conceivable that the preformed pair scenario of Fig. 1 is 
applicable. On the other hand we expect the PG to be 
present at large temperatures as a non-perturbative local 
lattice effect also in the two-dimensional lattice model. 
Certainly, other effects like strong phase fluctuations and 
small momentum pairing fluctuations, not contained in 
our calculations, can lead to an extension of the regimes 
where PG behavior occurs. 

A word of caution is in order when discussing the large 
temperatures addressed in this paper. Here we dealt with 
a strict one-band model where the kinetic energy is lim¬ 
ited by the bandwidth. In most real systems very high 
temperature would activate higher bands, and in solid 
state systems, it can lead to the melting of the crystal 
structure; such effects are obviously not allowed in our 


Appendix 

1. T-matrix approximation 


A popular approximation for the self-energy is the so- 
called T-matrix approximation, which corresponds essen¬ 
tially to summing the scat tering processes in the particle- 
particle channel. One has) 26 32 

£« =TUJ2 e iUmV G(q, iLO m ), (16) 

m,q 


or equivalently, 




q 


duj p(q , oj)np(uj), 


(17) 


and 


^k( iu> n) = T^2e lUriri T(q,iu> m )G(q - k,iu m - iu n ), 

(18) 


m,q 


with tj —y 0. Here we defined 

U K (q, icj m ) 


r(g, icj m ) — 


(19) 


1 - UK(q,iUm)’ 
with the particle-particle propagator 
K(q,iu m ) = -T^G(q-k,iu m -iu n )G{k,iuj n ). (20) 


n,q 


The self-energy is E fe (iw ra ) = + E 1(iu} n )- 

























In the local approximation, the expression simplify. 
We find the following result after analytic continuation, 


12 


= U Jdu (21) 


and 

S T (w) = IduJi lduj 2 [n B (^i) + n F (u> 2 )]- 

J J ^ U)2 

( 22 ) 

We have 

K(iuj m ) = -T y G(iuj m - iu) n )G(itv n ), (23) 


and pY = — ^Imr(u; + ). Introducing spectral functions 
we can also write, 


K(lo + ) = /dwi / du>2 [n F (o;i) —n F (—o; 2 )], 

J J ^ — CJi — CJ2 

(24) 


r(icj m ) 


U 2 K{iuj m ) 

1 - UKiiUm)' 


(25) 


The T-matrix calculations can be done non-self- 
consistently (Tnsc) and self-consistently (Tsc). 


2. Comparison of NRG-DMFT with IPT and 
T-matrix 


We start with a comparison of the DMFT results ob¬ 
tained using NRG calculations for the effective impurity 
model with DMFT calculations using second order per¬ 
turbation theory, usually termed iterated perturbation 
theory (IPT). IPT gives qualitatively reliable results in 
the half filled Hubbard modelP^ Since IPT does not re¬ 
quire a prescription of broadening discrete excitations, 
this comparison helps to validate the finite temperature 
broadening procedure described in Sec. II. We focus on 
results at half filling in this section. In Fig. [l3]we show 
a comparison of the imaginary part of the self-energy, 
Imy(cj), and the integrated spectral function, p(cj), for 
U/W = 0.6 (left) and U/W = 1 (right) and different 
temperatures. 

Overall the agreement is good with minor deviations in 
the tails. There is a particularly visible difference for 
U/W = 1, where the IPT result for Imy(cj) shows a 
somewhat stronger peak. This leads to a more pro¬ 
nounced PG in p{uj). We conclude that the DMFT-NRG 
results at high temperatures have the qualitative correct 
form and the PG remains there. 

We also provide a comparison of the DMFT-NRG re¬ 
sults with T-matrix calculations. In particular, we use 
the Eq. (21) and following, and include self-consistent 



FIG. 13: (Color online) Comparison of DMFT-NRG (full 
lines) and IPT (dashed lines) results for U/W — 0.6 (left) 
and U/W — 1 (right) for ImT(cj) and p{uS). 


(Tsc) and non-selfconsistent (Tnsc) results. Note that 
the T-matrix calculations are only sensible as long as 
1 — UReK(uj) does not become zero, which is particu¬ 
larly important for the non-selfconsistent case. At weak 
coupling ( U/W = 0.2, not shown) one can find reason¬ 
able agreement of T-matrix calculations with the DMFT- 
NRG and all calculations give no PG behavior. However, 
in this situation also second order perturbation theory 
gives satisfactory agreement. 



FIG. 14: (Color online) Comparison of DMFT-NRG result 
with self-consistent (Tsc) and non-selfconsistent (Tnsc) T- 
matrix calculations for ImE(cj) and p(uS) for U/W — 0.4 and 
T/W = 0.1. 


As seen in Fig. 14 for U/W = 0.4 and T/W = 0.1, 
Tsc and DMFT still show reasonable agreement, whereas 
Tnsc calculations can lead to a strong overestimate for 
ImE(u;). This can lead to a PG feature in p(o;), even 
though calculations with the DMFT-NRG give no PG 
behavior. 

For intermediate coupling, U/W = 0.6, and T/W = 
0.2, we show a further comparison in Fig. [l5j 
In this case both T-matrix calculations give unreliable re¬ 
sults. The self-energy of the self-consistent version is too 
small and p{uj) shows no PG. The non-selfconsistent cal¬ 
culation shows a PG but its magnitude is largely overes- 
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mensional Hubbard model at half filling. 


3. Second order self-energy and phase space factor 


The result for the second order retarded self-energy 
reads p 1 


FIG. 15: (Color online) Comparison of DMFT-NRG result 
with self-consistent (Tsc) and non-selfconsistent (Tnsc) T- 
matrix calculations for ImE(cj) and p(uj) for U/W = 0.6 and 
T/W = 0.2. 

timated. For larger interactions, for instance, U/W = 1, 
the deviations get worse. We therefore conclude that 
T-matrix calculations - both self-consistent and non- 
selfconsistent - within the local approximation do not 
give reliable results for the PG physics of the three di- 


V(w, k) = U 2 [de . (26) 

J u + irj-e 

The imaginary part of the retarded self-energy is then 
given by 

Im sr h (u) = -7rC/ 2 F r (w, k), (27) 

where F r (s, k) = /i(e, fe) + / 2 (f. k), with the phase space 
factors, 


/i(e,fe)= ^2 s (£k 2 +tk 3 -tk 1 -e)S(k + k 1 ,k 2 + k 3 )n kl (l-n k2 )(l-n k3 ). (28) 

k 1 ,k 2 ,k 3 

and, 

/ 2 (e,fc)= s (£k 2 +£k 3 8(k + k 1 ,k 2 + k 3 )(l-n kl )n k2 n k3 . (29) 

k i,k 2 ,k 3 


The expressions can be simplified in the limit of large dimensions. The momentum integrations can be replaced by 
integrals over the density of states, momentum conservation is implicit so we can omit the corresponding ^-function 
and the /c-dependence disappears, 


/i(e) = Jde 1 Jds 2 Jde 3 M £ i)M £ 2 )po( £ 3 )<K e 2 + e 3 -e 1 -e- p)n F (si - p)n F {-e 2 + p)n F (-e 3 + n). 


(30) 


We can do the integration over the ^-function, 


/i( £ ) 


= Jde 2 jde 3 p 0 (e 2 + e 3 


- e - ii)p 0 (£ 2 )pa{e 3 )n F (£2 + e 3 


e - p)n F (-e 2 + p)n F (-e 3 + p), 


(31) 


and similarly for / 2 (e). In the particle-hole symmetric 
case we have, 


/2(£)=/r(-£), (32) 

It is then sufficient to evaluate fi(s) and we can write, 

F*(e,k)=F r (e) = f 1 (e) + f 1 (-e). (33) 

This can be evaluated as a double integral for a given 
temperature and po(e). Assuming that po(e) is only fi¬ 
nite in an interval (—D,D) we can analyze the double 
integration as being determined by certain region in the 


£3 — £2 plane. At T = 0 a geometric analysis of the inte¬ 
gration region shows, fi(s) ~ £ 2 , which gives the typical 
Fermi liquid behavior, Eq. (10), at low temperature. In 
the opposite limit, T —)> 00, a similar analysis shows that 
F r (s) is maximal at e = 0 and it decays for small 5 as 


-£ , which yields the NFL form Eq. (11). One can es¬ 


timate the crossover temperature Tfl by studying when 
then coefficient of the e 2 changes sign. Depending on 
the density of states and the approximations made one 
finds a result of the order of a fraction of the bandwidth, 
consistent with the result in Fig. [2] for small U. 









14 


robert .peters@riken .j p 
^ jbauer@physics.harvard.edu 

1 T. Timusk and B. Statt, Reports on Progress in Physics 
62, 61 (1999). 

2 A. A. Kordyuk, ArXiv e-prints (2015), 1501.04154. 

3 S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. Nayak, 
Phys. Rev. B 63, 094503 (2001). 

4 D. J. Scalapino, Rev. Mod. Phys. 84, 1383 (2012). 

5 V. Emery and S. Kivelson, Nature 374, 434 (1995). 

6 V. Mishra, U. Chatterjee, J. Campuzano, and M. Norman, 
Nature Physics (2014). 

7 K. Efetov, H. Meier, and C. Pepin, Nature Physics 9, 442 
(2013). 

8 J. Bardeen, L. Cooper, and J. Schrieffer, Phys. Rev. 108, 
1175 (1957). 

9 D. M. Eagles, Phys. Rev. 186, 456 (1969). 

10 A. J. Leggett, in Modern Trends in the Theory of Con¬ 
densed Matter , edited by A. Pekalski and R. Przystawa 
(Springer, Berlin, 1980). 

11 P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1985). 

12 M. Randeria, in Bose-Einstein Condensation , edited by 
A. Griffin, D. Snoke, and S. Strinagari (Cambridge Uni¬ 
versity Press, Cambridge, 1995). 

13 M. Greiner, C. Regal, and D. Jin, Nature 426, 537 (2003). 

14 M. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Rau- 
pach, A. J. Kerman, and W. Ketterle, Phys. Rev. Lett. 92, 
120403 (2004). 

15 M. Zwierlein, J. Abo-Shaeer, A. Shirotzek, C. H. Schunck, 
and W. Ketterle, Nature 435, 1047 (2005). 

16 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 

17 J. K. Chin, D. E. Miller, Y. Liu, C. Stan, W. Setiawan, 
C. Sanner, K. Xu, and W. Ketterle, Nature 443, 961 
(2006). 

18 M. Feld, B. Frohlich, E. Vogt, M. Koschorreck, and 
M. Kohl, Nature 480, 75 (2011). 

19 J. Gaebler, J. Stewart, T. Drake, D. Jin, A. Perali, P. Pieri, 
and G. Strinati, Nature Physics 6, 569 (2010). 

20 M. Randeria, Nature Physics 6, 561 (2010). 

21 Q. Chen, K. Levin, and J. Stajic, J. Low Temp. Phys. 32, 
406 (2006). 

22 T. Eckl, D. J. Scalapino, E. Arrigoni, and W. Hanke, Phys. 
Rev. B 66, 140510 (2002). 

23 A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

24 Q. Chen, J. Stajic, S. Tan, and K. Levin, Physics Reports 
412, 1 (2005), ISSN 0370-1573. 

25 C.-C. Chien, H. Guo, Y. He, and K. Levin, Phys. Rev. A 

81, 023622 (2010). 

26 M. Keller, W. Metzner, and U. Schollwock, Phys. Rev. B 
60, 3499 (1999). 

27 D. Rohe and W. Metzner, Phys. Rev. B 63, 224509 (2001). 

28 F. Marsiglio, P. Pieri, A. Perali, F. Palestini, and G. C. 
Strinati, Phys. Rev. B 91, 054509 (2015). 

29 A. Perali, P. Pieri, G. C. Strinati, and C. Castellani, Phys. 
Rev. B 66, 024510 (2002). 

30 S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 
80, 033613 (2009). 

31 R. Watanabe, S. Tsuchiya, and Y. Ohashi, Phys. Rev. A 

82, 043630 (2010). 


R. Haussmann, Z. Phys. B 91, 291 (1992). 

R. Haussmann, M. Punk, and W. Zwerger, Phys. Rev. A 
80, 063612 (2009). 

M. Bauer, M. M. Parish, and T. Enss, Phys. Rev. Lett. 
112, 135302 (2014). 

P. Magierski, G. Wlazlowski, A. Bulgac, and J. E. Drut, 
Phys. Rev. Lett. 103, 210403 (2009). 

H. Hu, X.-J. Liu, P. D. Drummond, and H. Dong, Phys. 
Rev. Lett. 104, 240407 (2010). 

G. Wlazlowski, P. Magierski, J. E. Drut, A. Bulgac, and 
K. J. Roche, Phys. Rev. Lett. 110, 090401 (2013). 

M. Randeria, N. Trivedi, A. Moreo, and R. T. Scalettar, 
Phys. Rev. Lett. 69, 2001 (1992). 

N. Trivedi and M. Randeria, Phys. Rev. Lett. 75, 312 
(1995). 

Y. Vilk and A.-M. Tremblay, Journal de Physique I 7, 1309 
(1997). 

S. Moukouri, S. Allen, F. Lemay, B. Kyung, D. Poulin, 
Y. M. Vilk, and A.-M. S. Tremblay, Phys. Rev. B 61, 7887 
( 2000 ). 

A. Koga and P. Werner, Phys. Rev. A 84, 023638 (2011). 

A. Koga and P. Werner, Modern Physics Letters B 25, 973 

( 2011 ). 

M. Randeria and E. Taylor, Annual Review of Condensed 
Matter Physics 5, 209 (2014). 

B. Kyung, S. Allen, and A.-M. S. Tremblay, Phys. Rev. B 
64, 075116 (2001). 

M. Keller, W. Metzner, and U. Schollwock, Phys. Rev. 
Lett. 86, 4612 (2001). 

M. Capone, C. Castellani, and M. Grilli, Phys. Rev. Lett. 
88, 126403 (2002). 

B. Kyung, A. Georges, and A.-M. S. Tremblay, Phys. Rev. 
B 74, 024501 (2006). 

N. Kuleeva, E. Kuchinskii, and M. Sadovskii, Journal of 
Experimental and Theoretical Physics 119, 264 (2014). 

A. Garg, H. R. Krishnamurthy, and M. Randeria, Phys. 
Rev. B 72, 024517 (2005). 

A. Toschi, P. Barone, M. Capone, and C. Castellani, New 
J. Phys. 7, 7 (2005). 

A. Toschi, M. Capone, and C. Castellani, Phys. Rev. B 72, 
235118 (2005). 

J. Bauer and A. C. Hewson, Europhys. Lett. 85, 27001 
(2009). 

J. Bauer, A. C. Hewson, and N. Dupuis, Phys. Rev. B 79, 
214518 (2009). 

R. Micnas, J. Ranninger, and S.Robaszkiewicz, Rev. Mod. 
Phys. 62, 113 (1990). 

R. Bulla, T. Costi, and T. Pruschke, Rev. Mod. Phys. 80, 
395 (2008). 

R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev. B 
74, 245114 (2006). 

A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 
076402 (2007). 

J. Bauer, A. Oguri, and A. Hewson, J. Phys.: Cond. Mat. 
19, 486211 (2007). 

T. Hecht, A. Weichselbaum, J. von Delft, and R. Bulla, J. 
Phys.: Cond. Mat. 20, 275213 (2008). 

R. Bulla, T. A. Costi, and D. Vollhardt, Phys. Rev. B 64, 
045103 (2001). 

A. Georges, Annalen der Physik 523, 672 (2011), ISSN 
1521-3889. 


32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 



15 


63 We would like to thank N. Dupuis for pointing this out to 
us. 

64 L. De Leo, J.-S. Bernier, C. Kollath, A. Georges, and V. W. 
Scarola, Phys. Rev. A 83, 023606 (2011). 

65 E. Khatami, E. Perepelitsky, M. Rigol, and B. S. Shastry, 
Phys. Rev. E 89, 063301 (2014). 

66 We use this as a general condition also at finite tempera¬ 
ture, even though the Fermi surface is strictly only defined 
at T = 0. 


67 A. Kampf and J. R. Schrieffer, Phys. Rev. B 41, 6399 
(1990). 

68 A. P. Kampf and J. R. Schrieffer, Phys. Rev. B 42, 7967 
(1990). 

69 J. Bauer, M. Babadi, to be published (2015). 

70 J. Bauer and A. C. Hewson, Eur. Phys. J. B 57, 235 (2007). 

71 H. Schweitzer and G. Czycholl, Z. Phys. B 83, 93 (1991). 



